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Abstract. Hyperaccreting disks around neutron stars or magnetars cooled via neutrino emission can 
be the potential central engine of GRBs. The neutron-star disk can cool more efficiently, produce 
much higher neutrino luminosity and neutrino annihilation luminosity than its black hole counter- 
part with the same accretion rate. The neutron star surface boundary layer could increase the annihi- 
lation luminosity as well. An ultra relativistic jet via neutrino annihilation can be produced along the 
stellar poles. Moreover, we investigate the effects of strong fields on the disks around magnetars. In 
general, stronger fields give higher disk densities, pressures, temperatures and neutrino luminosity; 
the neutrino annihilation mechanism and the magnetically-driven pulsar wind which extracts the 
stellar rotational energy can work together to generate and feed an even stronger ultra-relativistic jet 
along the stellar magnetic poles. 
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INTRODUCTION 

The hyperaccreting disk surrounding a stellar-mass black hole possibly formed by the 
merger of a compact object binary or the collapse of a massive star has been argued to be 
a candidate for central engine of gamma-ray bursts (GRBs; e.g., Popham et al. 1999). On 
the other hand, both observational and theoretical evidences show that newborn neutron 
stars or magnetars rather than black holes can form in the GRB central engines (Usov 



1992; Dai & Lu 1998). However, all the previous magnetized neutron star or magnetar 
models ignore the accretion process which occurs onto a protoneutron star for the first 
several seconds. In the collapsar scenario, if the rotational core collapse can lead to the 
formation of a neutron star or a magnetar, it is possible that the prompt accretion or 
fallback process can make a hyperaccreting disk around the young formed star. Also a 
debris disk around a massive neutron star can be formed in compact binary mergers. 
Therefore we investigate the hyperaccreting neutron star or magnetar system, and show 
that it can also proposed as another possible central engine of GRBs (Zhang & Dai 
2008, 2009, 2010). For simplicity we consider the vertically integrated disk using the 
a -prescription. To our knowledge, our work is the first to study the hyperaccreting disks 
around neutron stars possibly related to GRBs. 



NEUTRINO-COOLED DISKS AROUND NEUTRON STARS 



In the hyperaccreting disk the Alfven radius of the central star can be estimated as 
r A ~ 0.201MZ^ 7 M ] ~y 7 jU3^ 7 km, where M = 0.01M- 2 M Q s _1 is the accretion rate, 
M = \AM\aMq is the central star mass, /i = /i 30 10 30 G cm 3 is the central magnetic 
flux. If the surface field B is less than B < 5 crit = 0.89 x 10 15 M^M^r^ 4 G with 

r = 5 6l0 6 cm being the star radius, the accretion flow will continue to be confined in 
the disk plane without co-rotating with the compact object or getting funneled onto the 
magnetar poles. 

In the weak field case, the quasi-steady disk around a neutron star can be approxi- 
mately divided into two regions — inner and outer disks, depending on the energy trans- 
fer and emission in the disk. For the outer disk, the heating energy rate Q + is mainly 
due to local dissipation (Q + = Q^ is ), and the structure of the outer disk is very simi- 
lar to the black hole disk. On the other hand, the heating energy in the inner disk in- 
cludes both the energy generated by itself and the energy advected from the outer region 
(Q + = Qvi s + QjjdyX so the inner disk has to be dense with a high pressure. We take 
Q + = Q and the entropy-conservation self- similar condition to describe the inner disk. 
The left panel of Figure 1 shows the size of the inner disk, which is determined by the 
global energy equation of the inner disk. Moreover, we also calculate the size of a low-oc 
inner disk, which is smaller for a low accretion rate (< O.1M s _1 ) compared to a high- 
er disk (a ~ 0.1), and increases dramatically with increasing accretion rate. Also, the 
inner disk structure would not exist if the disk is unstable to drive a large scale thermal 
non-relativistic wind when the accretion rate > O.5M s . 

We find that, due to the inner disk structure as well as the inner surface boundary of the 
compact star, the disk has a denser, hotter inner region with higher pressure compared 
to its black hold counterpart. Also, the entire disk can cool more efficiently via neutrino 
emission. The right panel of Figure 1 shows the total neutrino emission luminosity of 
the entire disk around a neutron star as a function of accretion rate, and we compare 
it with the neutrino luminosity from a black-hole disk. We see that the difference in 
neutrino luminosity between the neutron-star and black-hole cases is a strong function 
of the accretion rate. When the accretion rate is low, the total neutrino luminosity of the 
black-hole disk L v .bh is much smaller than that of the neutron-star disk L v .ns, but L v .bh 
and L V; ns are similar for a moderate accretion rate. 

NEUTRINO ANNIHILATION 

Hyperaccreting black hole disks can convert some fraction of the net accretion energy 
into the energy of a relativistic jet by two general mechanisms: neutrino annihilation 
and magnetohydrodynamical (MHD) effects such as the Blandford-Znajek mechanism. 
However, for hyperaccretion disks surrounding neutron stars, the energy conversion 
mechanism is mainly due to the neutrino annihilation for the magnetic fields < 10 15 
G. As the neutron star disk has a brighter neutrino luminosity compared with the 
black hole disk, the neutrino annihilation efficiency of the neutron star disk should be 
higher than that of the black hole disk. Moreover, the surface boundary of the neutron 
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FIGURE 1. (a) Left panel: the radius rj (10 6 cm) between the inner and the outer disks with the center 
neutron star M = 1 .4M Q (solid line), and M = 2.OM (dashed line), (b) Right panel: Neutrino luminosity 
from the disk with M = 1 AM & . The solid line corresponds to the maximum energy release rate of the disk 
around a neutron star, the dashed line to the neutrino luminosity from the inner disk, the dotted line to the 
neutrino luminosity from the outer disk, and the the dash-dotted line to the neutrino luminosity from a 
black hole disk. 
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FIGURE 2. (a) Left panel: Neutrino annihilation luminosity L V v (thick lines) and total neutrino emission 
luminosity L v (thin lines) as functions of accretion rate. The solid lines correspond to the neutron star disk 
with the boundary layer emission efficiency TJ S = 0, the dashed lines to r\ s = 0.5, and the dotted lines to 
the black hole disk, (b) Right panel: The neutrino annihilation rate l V y along the stellar m agnetic pole as a 
function of height for accretion rate M — 0.5M Q s _1 in strong field case, where sin0 ~ \j r^jr^ measures 
the strength of funnel accretion, smaller sin0 corresponds to stronger field and funnel accretion. 



star, which carries away gravitational-binding energy by neutrino emission makes the 
neutrino annihilation luminosity of the neutron star disk be even higher. We follow the 
approximate method used by Popham et al. (1999) to calculate the annihilation rate, and 
we introduce the efficiency factor r} s to measure the energy emitting from the stellar 
surface as L s ~ r\ s GMM/ (4r*) with r* being the neutron star radius. 

The left panel of Figure 2 shows the total neutrino annihilation luminosity and the 
emission luminosity of a neutron-star disk with different surface boundary layer condi- 
tions (T] s =0 and 0.5). The annihilation luminosity of the neutron- star disk are brighter 
than those of a black hole disk with the same mass and accretion rate. If we study the 
neutrino annihilation from the entire disk without surface boundary emission (t| v =0), 
the difference between L V v : ns and L V v,fi# is more significant for a low accretion rate 



than for a high accretion rate. On the other hand, neutrino emission from the neutron 
star surface boundary layer (r]. v =0.5) makes the annihilation luminosity be more than 
one order of magnitude higher than that without boundary emission (r] s = 0). Therefore, 
a lower-spin neutron star (larger rj s ) with hyperaccreting disk around it could have an 
obviously higher annihilation efficiency than that of a higher-spin neutron star. 

EFFECTS OF STRONG MAGNETIC FIELDS 

Strong fields > 10 15 G from the center magnetar can play a significant role in affecting 
the disk properties and even changing the accretion process (Zhang & Dai 2010). We 
consider the magnetar field has a dipolar vertical component B z . The differential rotation 
between the disk and the magnetar will generate a toroidal field component and a 
relatively weak radial component B r . The generated field can have an open or closed 
configuration, depending on the disk's viscous turbulence, magnetic diffusivity and disk 
angular velocity. For the strong field environment, the quantum effects (Landau levels) 
and large scale field coupling play two competitive roles in changing the disk properties, 
the former to decrease the pressure, density and neutrino luminosity with increasing field 
strength, while the latter to increase them. However, in most cases the field coupling is 
more significant than the microphysical quantum effect. Note that in the hyperaccreting 
disks, the funnel accretion can only be important for extremely strong fields. The funnel 
flow will cover a ring-like belt of "hot spot" around the magnetar surface and emit 
thermal neutrinos. Because of the more concentrated emission from the "hot spot", the 
funnel accretion can accumulate even more powerful neutrino annihilation luminosity 
than the weak field case (see right panel of Figure 2). 

The neutrino annihilation process both from the magnetar surface and from the disk 
plane will be higher than that without fields. Moreover, the neutrino annihilation mecha- 
nism and the magnetic activity from the stellar surface (i.e., the pulsar wind mechanism) 
can work together to generate and feed an ultra-relativistic jet along the stellar magnetic 
poles. If the stellar spin period is sufficiently short (e.g., ~4ms for the field ~ 10 16 
G and M = 0.1M© s _1 ), the jet from the magnetar will be magnetically-dominated and 
mainly feeded by extraction the stellar rotational energy. If the magnetar spin period is 
long, the jet is thermally-driven and feeded by the annihilation process. In this case, an 
energetic ultrarelativistic jet via neutrino annihilation can be produced above the stellar 
polar region if the disk accretion rate and the stellar surface luminosity are sufficiently 
high. In the intermediate case, on the other hand, the relativistic jet can be launched by 
the pulsar-wind-like process and neutrino annihilation together. 
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